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^ ' Abstract. We estimate the total mass distribution of the galaxy cluster CL0024+1654 from the measured source 

depletion due to lens magnification in the 7? band. Within a radius of 0.54 h ^ Mpc, a total projected mass of 
(8.1 ±3.2) X lO"/i-^M0 (EdS) is measured. The la error here includes shot noise, source clustering, uncertainty 
in background count normalisation and contamination from cluster and foreground galaxies. This corresponds 
to a mass-to-light ratio of M/Ls = 470 ± 180. We compute the luminosity function of CL0024-I-1654 in order 
to estimate contamination of the background source counts from cluster galaxies. Three different magnification- 
based reconstruction methods are employed: 1) An estimator method using a local calculation of lens shear; 2) A 
non-local, self-consistent method applicable to axi-symmetric mass distributions; 3) A non-local, self-consistent 
method for derivation of 2D mass maps. We have modified the standard single power-law slope number count 
theory to incorporate a break and applied this to our observations. Fitting analytical magnification profiles of 
Q ' different cluster models to the observed number counts, we find that CL0024-I-1654 is best described either 

^ '. by a NEW model with scale radius Tb = 334 ± 191 kpc and normalisation Ks = 0.23 ± 0.08 or a power- 

C/3 . law profile with slope ^ = 0.61 ± 0.11, central surface mass density kq = 1.52 ± 0.20 and assuming a core 

' radius of rcom ~ 35 /i"'^ kpc. The NEW model predicts that the cumulative projected mass contained within a 

radius R scales as M(< R) = 2.9 x 10^''(i?/l')^'^"° '^ ''^^'^^^''/i'^Mq . Einally, we have exploited the fact that fiux 
magnification effectively enables us to probe deeper than the physical limiting magnitude of our observations in 
searching for a change of slope in the U band number counts. We rule out both a total flattening of the counts 
with a break up to Uab < 26.6 and a change of slope, reported by some studies, from dlogA^/dm = 0.4 0.15 
up to Uab < 26.4 with 95% confldence. 

Key words. Cravitational lensing - Galaxies: clusters: individual: CL0024-f 1654 - dark matter 



1. Introduction 

The lensing cluster CL0024-fl654 ranks as one of the 
most highly studied clusters to date. Lying at a redshift 
of z = 0.39, early measurements of the cluster's velocity 
dispersion of cr ~ 1300± 100km s^^ (Dressier, Schneider & 
Gunn 1985) suggested a formidable mass. The discovery of 
a large gravitationally lensed arc from a blue background 
galaxy by Koo (1988) has since provoked a range of studies 
to estimate the cluster's mass based on its lensing prop- 
erties. 

The first lens inversion of CL0024-I-1654, by Kassiola, 
Kovner & Fort (1992), noted a violation of the 'length 
theorem' (Kovner 1990) that the length of the middle seg- 
ment of the arc should equal the sum of the other two. 
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The authors demonstrated that a concentration of large 
cluster galaxies near the arc centre is necessary to cause 
this by perturbing the cluster cusp and were thus able 
to constrain the potential of the cluster and the perturb- 
ing galaxies. A later analysis by Bonnet, Mellier & Fort 
(1994) constrained the cluster's mass profile more tightly 
with the first measurement of weak shear out to a ra- 
dius of 1.5 h-^ Mpc. WaUington, Kochanek & Koo (1995) 
confirmed the perturbing galaxy hypothesis of Kassiola, 
Kovner & Fort by fitting a smooth elliptical cluster po- 
tential with two superimposed galaxy potentials near 
the centre of the arc. 

By parameterising the source and lens models and fit- 
ting to six images of the lensed source galaxy, Kochanski, 
Deir Antonio & Tyson (1996) again showed that the 
mass profile of CL0024-f 1654 is consistent with a smooth 
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isothermal distribution. Furthermore, they found that the 
cluster mass traces light fairly well out to a radius of 
0.5/i~^Mpc and were able to rule out the existence of 
any significant substructure larger than 15 h^^ kpc in the 
central region. Using HST images of the cluster, Tyson, 
Kochanski & DeU' Antonio (1998, TKD hereafter) isolated 
eight well-resolved images of the blue background galaxy 
to construct a high resolution mass map of the cluster. 
Their parametric inversion concluded that more than 98% 
of the mass concentration excluding that contributed from 
discrete galaxies was represented by a smooth distribution 
of mass with a shallower profile than isothermal. 

The most recent lensing analysis of CL0024+1654 to 
date is that by Broadhurst et al. (2000, BOO hereafter) 
who provide the first measurement of the redshift of the 
blue background galaxy dX z — 1.675. This redshift breaks 
the mass-redshift degeneracy present in all mass estimates 
of the cluster thus far. Their fit of NFW profiles (Navarro, 
Frenk & White, 1997) to the eight brightest cluster mem- 
bers is found to be an adequate model to explain the po- 
sitions of the five main lensed images. This, they suggest, 
highlights the possibility that sub-structure has not been 
erased in the cluster. 

The first of only two measurements of the cluster's 
lens magnification to currently exist was that investi- 
gated by Fort, Melher & Dantel-Fort (1997). This was 
the first detection of depletion in background galaxy num- 
ber counts due to geometrical magnification as predicted 
by Broadhurst, Taylor & Peacock (1995, BTP hereafter). 
Rather than reconstruct cluster mass, this work concen- 
trated on characterising the radial distribution of critical 
lines to infer the redshift range of the background popu- 
lations in B and /. Using this data, van Kampen (1998) 
produced an estimate of the mass of CL0024-I-1654 from 
lens magnification. 

The only other magnification analysis of the cluster 
published so far is that of Rognvaldsson et al. (2001; ROl 
hereafter) using R and U band observations. Their choice 
to observe in U was inspired by the findings of Williams 
et al. (1996) which suggested a fiattening in the U band 
number count slope at faint magnitudes. Such a break is 
reported to occur at Uab — 25.5 — 26 with a change of 
slope of d log TV/dm ~ 0.4 0.15. Given suitably deep 
U band imaging, this should therefore manifest itself as 
a depletion in the number density of galaxies observed in 
the presence of lens magnification (see Section 3.2). ROl 
claimed to have detected depletion in U implying that a 
break in the slope must be present. 

Further investigation has since highlighted concern re- 
garding the reliability of faint objects extracted by ROl 
from the R and U band observations. While this is not a 
large concern for R where the depletion signal is strong, 
the claim of detection of depletion in U and hence the 
reported break in the number count slope is strongly af- 
fected. A re-examination of these findings is therefore nec- 
essary. 

The motivation driving the paper presented here is sev- 
eralfold. Firstly, there is the need to re-evaluate the exist- 



ing results of ROl. We take the observations from ROl 
but create new, more reliable R and U band object cat- 
alogues. Using these, we re-calculate the depletion signal 
in both bands and use this to re-fit the isothermal lens 
model of ROl. Secondly, we wish to extend the analysis 
of ROl: 1) We fit a power-law profile and a NFW model 
to the depletion profiles. 2) We transform the measured 
depletion into mass estimates using three recently devel- 
oped magnification reconstruction methods; the local es- 
timator method of van Kampen (1998) applied by Taylor 
et al. (1998, T98 hereafter), the non-local axi-symmetric 
solution of T98 for reconstruction of radial mass profiles 
and the non-local method of Dye & Taylor (1998) for de- 
termining 2D mass maps. 3) We investigate the relation- 
ship between mass, light and galaxy number density in 
the cluster. 4) We quantify contamination of our back- 
ground source sample by cluster and foreground objects 
and incorporate this into our analyses. 5) Finally, we apply 
two methods to search for a change in slope in the U band 
field galaxy number counts by exploiting the fact that lens 
magnification enables us to effectively see deeper than the 
physical magnitude limit imposed by the observations. 

The following section briefiy describes data acquisi- 
tion, reduction, object extraction and mask generation. 
Contamination of the source counts from cluster and fore- 
ground galaxies is estimated. Section ^ details the recon- 
struction methods used in this paper. In Section |4[ we 
consider the properties of the R band background galaxy 
population and fit cluster mass models to the measured 
number count profile. This is necessary for determination 
of the mass results presented in Section ^. In Section |[ 
we investigate the U band galaxy population and test for 
the existence of a break in the number counts. Finally, we 
discuss and summarise our findings in Section |^. 



2. The data 

Data acquisition and reduction are described in full detail 
in ROl. Here, we highlight the key deviations in the gen- 
eration of the object catalogues of this paper from those 
in ROl: 

— A Gaussian smoothing kernel with a FWHM equal to 
the seeing was used to extract objects (see Section 2.2). 
This compares with the narrower smoothing kernel of 
FWHM equal to slightly larger than half the seeing 
used by ROl (note the incorrect statement that their 
smoothing kernel is equal to the seeing). Matching the 
width of the smoothing kernel to the seeing allows op- 
timal extraction in the sense that spurious detections 
due to background noise are kept as small as possible 
while maintaining a sufficiently large detection success 
rate. This creates the most noticeable difference be- 
tween this dataset and that of ROl; our U and R cat- 
alogues contain fewer faint objects (our total paired 
object catalogue contains 29% fewer objects). 



Dye et al.: Lens magnification by CL0024+1654 in the U and R band 



3 



— Galactic extinction corrections have been applied to 
both the U and R band, brightening all objects (see 
Section p]). 

— A more thorough determination of the background 
noise has been obtained. The mean RMS background 
variation has been calculated within a circular aper- 
ture of diameter equal to the seeing FWHM. The re- 
sult of this is that the 'Sfi' detection limits calculated 
in ROl from Poission statistics are fainter than those in 
this paper by 1.3 mag in U and 1.7 mag in R including 
the galactic extinction correction. 

— A linear co-ordinate transformation has been applied 
to the U band object positions resulting in a more ac- 
curate mapping to the R band. This has improved the 
position coincidence matching between both bands. 

2.1. Acquisition and completeness 

CL0024-I-1624 was observed in the Cousins U and Cousins 
R bands with a 6' x 6' field of view using the Nordic Optical 
Telescope in August 1999. A total integration time of 37 
ksec and seeing 1.1" was obtained in the U band com- 
pared with 8.7 ksec and 1.0" in the R. Following ROl, we 
use AB magnitudes throughout this paper converting from 
the Cousins U and R magnitudes with an offset of -1-0.71 
mag and -1-0.20 mag respectively (Fukugita et al. 1995). In 
addition, we have corrected the data for galactic extinction 
amounting to —0.15 mag and —0.31 mag in the R and U 
band respectively (Schlegel et al. 1998). The limiting mag- 
nitudes corresponding to a signal-to-noise ratio of 3 within 
a seeing disk are C/ab = 25.7 mag and i?AB = 25.8 mag, 
with the noise level taken as the mean RMS of the back- 
ground. All photometry was performed using SExtractor 
(Bertin & Arnouts 1996; see Section |2!^) . Magnitudes were 
measured as either total or corrected isophotal depending 
on the proximity of neighbouring objects. This is man- 
aged automatically via SExtractor 's MAG_BEST magnitude 
definition. 

Completeness was estimated from the detectability of 
synthetic objects of varying brightness added to the im- 
ages. Further details can be found in ROl. In re-applying 
this process to the data of this paper, we find that the com- 
pleteness at the 3tT detection limit is 84% at Uab = 25.7 
mag and 81% at i?AB = 25.8 mag. 

2.2. Object extraction 

Objects were extracted from the final reduced images us- 
ing SExtractor. With a detection threshold of la above 
background and a Gaussian filtering kernel of FWHM 
equal to the seeing, catalogues of all objects with at least 
10 connecting pixels brighter than the threshold were gen- 
erated. A total of 1887 objects in the R band and 1122 
objects in the U band were detected. These totals include 
stars but exclude objects classified by SExtractor as hav- 
ing saturated pixels, being truncated or possessing cor- 
rupt isophotal data. By matching the 30 brightest stars 



between both bands, a linear co-ordinate transformation 
was calculated and applied to map the U band catalogue 
onto the R band. Objects were paired within a positional 
tolerance equal to the seeing and yielded a total of 875 
objects excluding stars. Those objects with star/galaxy 
classification indices larger than 0.95 were assumed to be 
stars and excluded from the analysis below. 

2.3. Object selection 

Separation of the background galaxies from the foreground 
and cluster galaxies must be achieved before lens magni- 
fication can be evaluated. Segregation of the cluster and 
foreground objects is important to enable estimation of 
their background sky obscuration which affects binned 
number counts. 
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Fig. 1. Colour-magnitude diagrams for all 875 matched ob- 
jects. Top and bottom shows selection of mask objects in the 
U and R band respectively by the criteria {U — R)ab > 3 and 
Uab < 24, Rab < 23.5 (dashed hnes). The 3a detection hmits 
Uab = 25.7 and -Rab = 25.8 are shown in both plots by the 
dot-dashed hnes. 

Figure ID shows the colour-magnitude plot of {U—R)ab 
versus the Uab a-nd i?AB magnitudes for the 875 matched 
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objects. By matching our sample with the redshift mea- 
sm'ements of the field of CL0024+1654 by Czoske et al. 
(2001b), we have found that a large fraction of cluster 
members can be immediately discarded by removing ob- 
jects with ([/ — R)ab > 3. The noticeable lack of objects 
seen in the vicinity of {U — R)ab — 3 reflects a mini- 
mum in the bi-modal [U ~R) distribution of cluster galax- 
ies identified in the Czoske et al. sample by the criterion 
0.388 < z < 0.405. Applying a further selection i?AB < 22 
to avoid incompleteness in both the Czoske et al. sample 
and our data, we find that 65% of identified cluster ob- 
jects lie at iU — R)ab > 3. In addition, we find that 12% 
of objects brighter than i?AB — 22 with {U — i?)AB > 3 
are foreground galaxies and only 2% of objects in this se- 
lection are background objects. In summary, the selection 
{U — R)ab > 3 very efficiently removes cluster galaxies, 
at least up to Rab = 22, discarding a very small fraction 
of background sources. 

In Figure H, we plot the number density of cluster 
galaxies defined as objects with colour {U — R)ab > 3 
or by their redshift for {U — R)ab < 3. Figure ^ shows 
how the number counts of galaxies meeting these crite- 
ria vary with magnitude in the U and R. In Figure ^ we 
plot the distribution of light from these selected cluster 
galaxies. The surface fiux density shown in this plot is 
scaled to the restframe B band for easy comparison of the 
mass-to-light ratio we calculate later with other authors. 
This scaling was carried out using the K and colour cor- 
rections presented in Fukugita et al. (1995). Both distri- 
butions are com par ed with the reconstructed cluster mass 
map in Section |5.l| . The fact that both distributions arc 
relatively concentrated around the known cluster centre 
indicates that we have reliable selection criteria for the 
(bright) cluster members. 
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Fig. 2. Number density (n=number/square arcmin) of cluster 
galaxies selected either by (?7 — 7?)ab > 3 or redshift. 



X (arcmin) 

Fig. 3. Distribution of light (surface flux density L con- 
verted to the restframe B band and expressed in units of 
IO^^LqMpc"'^ h^) from cluster galaxies selected either by {U — 
R)ab > 3 or redshift. Contours are separated by equal fiux 
intervals. 



Our foreground and cluster objects are therefore cho- 
sen as a combination of objects identified in the Czoske 
et al. sample with z < 0.405 as well as those satisfying 
{U — R)ab > 3 and an additional i?AB < 23.5 for the R 
band sample or an additional Uab < 24.0 for the U band 
sample. The choice of R and U limits here are chosen to be 
optimal in the sense that they select as many background 
objects as possible while preventing too high a degree of 
contamination from cluster and foreground objects (see 
Section ^ ). Extracted parameters such as size, elliptic- 
ity and orientation of the cluster and foreground objects 
are used to generate an obscuration mask as Section 2.5 
details. 

Rather than perform our analyses on the associated 
object catalogue (ie. containing only objects with paired 
U and R mags) which would miss faint objects detected 
in only either the U or R band, we take our background 
source list for each band from the individual complete 
U and R catalogues. To obtain our source list, we re- 
move objects from each complete catalogue identified as 
cluster members and foreground objects using the colour- 
magnitude and redshift criteria above. 

For our U band break analysis in Section ^, removal of 
contaminants causes a slight dilemma. Although we would 
like our U band sample to be contamination free so that 
any lensing signal present is maximised, we also wish to 
remain consistent with existing U band studies to pre- 
vent our search for the break from being biased in some 
way. It turns out that since the Czoske et al. redshift sur- 
vey extends to only relatively shallow U band magnitudes, 
we only identify and remove 5 foreground objects which 
fall within the colour-magnitude selection criteria for our 
background sample. We have therefore introduced only a 
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Fig. 4. Number counts as a function of magnitude for the U 
band (top) and R band (bottom). Plotted are counts of galaxies 
in the background sample (solid) and cluster galaxies (dashed) . 
Errors account for shot noise only. 



negligible inconsistency with other U band number count 
studies and yet have removed all we can in terms of known 
foreground objects. We know that there are more fore- 
ground objects in our background sample than we have 
been able to remove (see Section 2.4.1) so we are forced 
to allow for these as a contamination error in our lensing 
analysis. As far as removal of objects with {U — R)ab > 3 
is concerned, this simply removes the majority of surplus 
cluster galaxies which aren't present in published number 
count studies anyway. 

In the R band, this does not pose a problem. We do 
not have to ensure that our R band data is consistent 
with any other sample since our depletion is normalised to 
counts near the edge of the field rather than independent 
measurements. Removal of as many foreground and cluster 
objects as possible is therefore the ideal. As with the U 
band sample, we can not achieve this fully so must resort 
to treating the remainder as a source of error in the mass 
reconstructions. 

After removal of cluster and foreground objects from 
our full U and R catalogues, we find a total of 863 back- 
ground galaxies remaining in the U and 1367 in the R. 



Figure ^ plots the number counts of both these samples 
as a function of magnitude. Notice how the U band counts 
are significantly steeper than in R. All subsequent magnifi- 
cation analysis will be performed on these two catalogues. 

2.4. Contamination of background source samples 

We discussed in the previous section that ~ 65% of clus- 
ter galaxies lie at (U — R)ab > 3 for Rab < 22. For the 
R band sample, this does not provide a useful means of 
removing cluster galaxies as the background source selec- 
tion limit of Rab > 23.5 prevents inclusion of objects with 
{U — R)ab > 3 anyway due to the U band detection limit 
of Uab — 25.7. By including objects from the full R band 
catalogue, faint objects with {U~R)ab > 3 will inevitably 
fall into the background R sample. Additionally, the clus- 
ter galaxies which lie at (C/ — -R)ab < 3 will contaminate 
both U and R samples along with any foreground galaxies. 
This contamination must be quantified. 

2.4.1. Foreground galaxy contamination 

To estimate the foreground galaxy contamination ex- 
pected within our chosen magnitude ranges, we use the 
luminosity functions measured by the CN0C2 field galaxy 
rcdshift survey (Lin et al. 1999). The CN0C2 survey is 
ideal for our purposes being the largest intermediate red- 
shift survey to date with multicolour U BVRI photometry. 

The number of objects N which exist within the ap- 
parent magnitude range mi < m < m2 and the redshift 
range zi < z < Z2 can be estimated as 

f-M2(z) 



N = 



dV{z) 



Mi{z) 



<j){M) dM 



(1) 



where dV{z) is the comoving volume element, 4>{M) is the 
luminosity function and the absolute magnitudes Mi{z) & 
M2{z) correspond to the apparent magnitudes mi and 7712 
at a redshift z. The fraction of objects which lie closer than 
a redshift in a sample of galaxies observed within the 
magnitude range mi < m < m2 is therefore given by, 

N(rni < m < m2, z < Zf) 
N{mi < m < 1112) 

The denominator here is the total number of galaxies 
within mi < m < m2 integrated over all redshifts. 
Using this equation, we estimate the fraction of galax- 
ies within our U and R background samples lying at 
z < 0.405, allowing for an early-|-intermediate mixed 
galaxy K-correction (Fukugita et al. 1995) for the mag- 
nitudes in equation (|l|). 

We use the early-|-intermediate (Cousins) R and 
(Cousins) U luminosity functions from the CN0C2 sur- 
vey for calculation of N in equation ([^). These predict 
that the contamination due to foreground galaxies is 3% 
in the U band sample and 2% in the R band assuming 
an EdS cosmology. For the case = 0.3, A = 0.7, these 
fractions drop only very slightly. We find that this result 
is insensitive to the choice of K-correction mix. 
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2.4.2. Cluster contamination and the cluster 
luminosity function 

To determine the cluster contamination fraction, we need 
to be able to predict the number of cluster members ex- 
pected within the magnitude range spanned by the se- 
lection magnitude and the detection limit in both bands. 
This requires knowledge of the cluster luminosity function 
(CLF) for CL0024+1654. 

We fit our own CLF to the cluster counts determined 
by matching with the Czoske et al. cluster members (see 
Section 2.3) in both the U and R. We use the likelihood 
method of Sandage Tammann & Yahil (1979) and check 
the consistency of our results with several published CLFs. 
In the R band, we compare with three recent studies of 
the composite CLF (CCLF): Paohllo et al. (2001) who 
construct a CCLF from 39 Abell clusters, Piranomonte 
et al. (2000) whose CCLF is constructed from 86 AbeU 
clusters and Garilli et al. (1999) who use 65 Abell and X- 
ray selected clusters to calculate a CCLF. All three studies 
use observations in the Gunn r filter so we must apply a 
K-colour-correction for the absolute magnitude conversion 
-^i?AB ~^ ■ Since ~ 60 — 70% of galaxies within clusters 
are E/SO types (Dressier et al. 1997; see also Small et al. 
1997 for specific case of CL0024+1654), we calculate this 
correction using an elliptical galaxy spectrum taken from 
Kinney et al (1996). 

In the U band, publications on CLFs are rare and there 
are no known U band CCLFs to date. Applying K-colour- 
corrections to convert redder-band magnitudes to the U 
band when the galaxy type mix is not known accurately is 
not reliable. We therefore choose the most heavily studied 
cluster in the U band, the Coma cluster, and compare our 
fit with the CLF from the recent study by Beijersbergen 
et al. (2001). In a similar manner to the R band, we apply 
a correction to take absolute Uab magnitudes at z = 0.39 
to rest-frame absolute Cousins U magnitudes. 

Figure || plots the various CLFs/CCLFs, all of which 
are described by a Schechter function (Schechter 1976) 
with parameters given in Table 0. For the r band CCLF 
from each study, we have taken Schechter parameters cor- 
responding to the rich subsample of clusters for consis- 
tency with CL0024-I-1654. All luminosity functions are 
normalised to our cluster counts which are superimposed 
in both plots. Absolute magnitudes assume qq = ho = 0.5. 
For correct normalisation, we must ensure that the Czoske 
et al. matched sample for each band does not suffer 
from incompleteness. We therefore limit the matched U 
band counts to the magnitude range 21.8 < Uab < 24.0 
(-22.4 < Mu < -20.2) and the matched R band counts 
to 19.9 < Rab < 22.1 (-23.3 < Af,. < -21.1). 

Integrating over each of the CLFs/CCLFs, we are able 
to calculate the number of cluster galaxies expected. For 
the published luminosity functions presented here, this is 
a straight-forward integral. For our own fitted CLFs how- 
ever, we incorporate the uncertainty arising from their fit 





Fig. 5. Luminosity functions of CL0024-I-1654 galaxies. Top: 
U band Schechter LF from Beijersbergen et al. (2001) for the 
Coma cluster (dashes) and the best fit LF determined in this 
work (solid). Bottom: r band CCLFs for rich clusters from 
Paolillo et al. (2001, dots), Piranomonte et al. (2000, dashes), 
Garilli et al. (1999, dot-dashes) and the best fit from this work 
(solid). 



to obtain a more fair level of contamination. We calculate 
the number of cluster galaxies as 



Mb \ 

ric^ A{(j)*{M^,a) I dM0(M;M*,a) 



(3) 



M,.a 



where the angular brackets denote averaging over each 
realisation of M^, and a weighted by the corresponding 
probability from the CLF fit. ^4 is a normalisation con- 
stant calculated as the inverse of the sum of these proba- 
bilities. The absolute magnitudes Ma and Mi, correspond 
to the apparent magnitude ranges 23.5 < Rab < 25.5 for 
r or 24.0 < Uab < 25.5 for Uc- These ranges are chosen 
to extend slightly less deep than the 3tT detection limits, 
allowing comparison with a more complete background 
sample to give a more accurate estimation of contamina- 
tion. 
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CLF/CCLF 



Band M. 



C (%) 



Paolillo CCLF 


r 


-22.2 


-1.01 


15 


Piranomonte CCLF 


r 


-22.0 


-0.91 


13 


Garilli CCLF 


r 


-22.1 


-0.60 


5 


This work CLF 


r 


-21.9 


-0.70 


19 


Beijersbergen CLF 




-19.4 


-1.54 


8 


This work CLF 




-21.5 


-1.41 


9 



Table 1. Schcchtcr parameters of the CLF/CCLFs used 
for estimation of cluster contamination, C, expressed as 
the percentage of objects expected to be cluster members 
within 23.5 < Rab < 25.5 for r or 24.0 < Uab < 25.5 for 
Uc- Note that the contamination fractions calculated from 
the CLFs of this work are weighted averages (see text). 



Within the appropriate magnitude ranges given above, 
the last column of Table ^ lists the predicted number of 
cluster galaxies as a percentage of the total number of 
galaxies detected. In R, the estimated contamination from 
the published luminosity functions ranges between 5% and 
15%. The weighted average predicted contamination from 
our own R band CLF is higher at 19%. In the U band, 
the contamination levels are not as high since the known 
Czoske cluster galaxies have been discounted in the cal- 
culation. We find that the Coma CLF predicts that 8% 
of objects in our background sample are cluster objects, 
compared with our CLF which gives a weighted average 
of 9%. 

The total contamination predicted from both fore- 
ground and cluster galaxies therefore has a relatively wide 
spread in the R band of 7% to 21% compared to between 
11% and 12% in U. Unfortunately, none of the published 
CLFs/CCLFs presented here extend sufficiently deep to 
be able to reliably predict the number of cluster objects 
in CL0024-I-1654 at the faint end of our sample. In addi- 
tion, some CLF studies (eg. Driver et al. 1994; Wilson et 
al. 1997) indicate the possibility of steeper faint end slopes 
which would increase the fraction of cluster contaminants. 
In light of this, we therefore assume the generous contami- 
nation levels predicted from our fitted CLFs in the analysis 
which follows in this paper. 

Clearly, the contamination is dominated by the cluster 
members. This means that the variation in contamination 
across the field of view is governed by their density profile. 
Since this increases toward the centre of the cluster and 
since the background number counts fall off toward the 
centre where the magnification is strongest (see Section 
4.3) the observed number density of objects in the back- 
ground sample is affected most by contamination at small 
radii. This distribution will have a strong fall-off as one 
moves away from the cluster centre. Fortunately, this re- 
duces the impact of the cont amin ation on the cumulative 
mass measurements (Section ^.3| ) at large radii where the 
strongest mass contribution comes from. 

An approximate contamination profile can be obtained 
by assuming that the cluster mass distribution and hence 
the number density of galaxies is described sufficiently ac- 



curately (for this purpose at least) by an isothermal sphere 
(see Section 4.3). The radial dependence of the error on 
the surface mass density, k, can thus be written 



(r) 



2nr(/3- l)/i/3(r) 



(4) 



where n is the background number density in the absence 
of lensing, /i(r) is the magnification given by equation (24) 
and f3 is the number count slope (see Section 3.1). We 
have approximated the number density profile of the clus- 
ter galaxies as k/r where k is set by normalising to the 
contamination fractions given earlier in this section. 

In all reconstructions of k found later in this paper, we 
add the error given by equation (^) in quadrature to the 
sum of the other errors, including the error due to fore- 
ground galaxy contamination. The effect of contamination 
on the search for the U band break is discussed separately 
in Section ||. 

2.5. Obscuration masks 



To account for obscuration of the background sky by clus- 
ter members and foreground objects, we create a mask. 
Without giving consideration to obscuration, the ratio of 
observed objects to expected objects used later in our 
analysis would be biased to lower values. This is especially 
true of bins near the cluster centre where relatively heavy 
obscuration occurs due to large central cluster galaxies. 

Figure ^ shows the obscuration mask produced us- 
ing the SExtracted parameters of mask objects selected 
The top half shows the U band mask 
U band background source positions (black 



in Section 2.3 



(grey ellipses) 

crosses) and the position of the annular bins used in 
Section |4.3| and Section ^ to bin counts radially. The lower 
half of Figure ^ shows the corresponding R band mask and 
background sources with the grid used in Section |5.1| for 
2D binning. Note the horizontal phase of the bins with 
respect to the field of view to allow placement of bin (7, 5) 
directly over the cluster centre, aligning it with the area 
encompassed by the critical line. 

3. Mass reconstruction theory 

In this section, we outline the process used for the deter- 
mination of magnification from number counts. We then 
detail the three different methods used to transform from 
magnification to surface mass density. 

3.1. Lens magnification: Single power-law number 
counts 

A background population of galaxies whose integrated 
number counts follow the standard power-law n{< m) oc 
2Q0.4/3m ^^Yi be observed under a lens magnification fac- 
tor fi to have the number count n'(< m) given by (BTP, 
T98); 



(5) 
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Fig. 6. Top: U band mask (grey ellipses) with sour ce po sitions 
(black crosses) and annular bins used in Section 4.3 and |^. 
Bottom: Corresponding R band mask with grid used for 2D 
maps in Section Observed arcs lie along the outer heavy 
dashed circle in both plots. The inner dashed circle shows the 
critical line determined from the isothermal sphere model fit 
to the R band number counts (see Section 4.3). 



The quantity Sni accounts for perturbations in n due to 
non-linear clustering and we define A as the expected num- 
ber of sources in the absence of clustering. We model the 
fluctuation Sni with a lognormal distribution and combine 
this with the additional uncertainty due to shot noise to 
give a Poisson-lognormal distribution (eg. Coles & Jones 
1991; BTP). Unlike a Gaussian distribution, the lognormal 
distribution accounts for non-linear clustering of the back- 



ground probed by our small field of view and is positive- 
definite. 

The Poisson-lognormal distribution can be defined as 
a compound distribution formed from a Gaussian distri- 
bution G of a linearised field and a Poission distribution 
D as 



PLN[«|A(Ai)] = / d6D[n\Xe 



V2 



]G{S). 



(6) 



Here, S is the linear density fluctuation which relates to 
the non-linear fluctuation as 1 -I- 6ni = e^~" with a the 
linear clustering variance. Following the method of T98 
but applied to the R band correlation function of Hudon 
& Lilly (1996), we calculate the non-linear clustering vari- 
ance. Averaging the angular correlation function over a 
circular area of radius 6 gives. 



1.5 X 10" 



-1.8 



(0/1') 



-0.8 



(7) 



The linear clustering variance is calculated from the non- 
linear variance using cr^ = ln(l + cr^^j). 

Equation (^ gives the probability distribution for the 
lens magnification in a given bin containing n galaxies. 
This directly gives the most probable magnification for 
that bin with its associated error taken as the width of the 
distribution. The integral in equation (^) must be evalu- 
ated numerically. 

3.2. Lens magnification: Dual power-law number 
counts 



In Section 3.1, we discuss evidence which suggests the 



presence of a break in the U band number count slope 



at faint magnitudes. Section 4.1 also highfights the possi- 



bility of a break in the R band counts. To account for the 
effect this has on the magnification calculated from num- 
ber count depletion, we modify equation ^ to incorporate 
a second slope applicable beyond some break magnitude 
nib. 

Writing the observed differential number counts as a 
dual power-law, 



, , , alO"-^''!'" , m < mf, 



(8) 



where the normalisation coefficients a and h are con- 
strained by continuity, the integrated number counts up 
to some limiting magnitude mi > rrih for fi > 1 can be 
expressed as 



n'{<mi) = fi [ni{< mb)+ 

n2i< mi)n'^^ - n2{< rrib)] 



(9) 



Here, ni and n2 denote number counts integrated up to 
a limiting magnitude over a constant slope of f3i and P2 
respectively. The clustering term which features in equa- 
tion (^ has been omitted here for clarity but is used in 
our analysis later. 
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If the limiting magnitude of observations matches the 
break magnitude then 



"-2(< m) = ^"-i(< 

Equation (0) can thus be simphfied to 



n = n/i 



1 + f (/^ 

P2 



(m>i) 



(10) 



(11) 



where n = ni(< m;) here, is the unlensed surface number 
count density observed, for example, at the edge of the 
field of view. If (3i = (32 and hence there is no break, 
equation (|ri|) becomes equation (||). For the case when 
/i < 1, equation (^) must be reverted to. 

3.3. Local mass estimator 

The simplest means of arriving at a mass estimate given 
a measurement of magnification is by assuming a local 
relationship between the convergence and shear. Such a 
relationship can be derived once a model is chosen to de- 
scribe the mass distribution. Since magnification depends 
on the convergence, k, and shear, 7, as 



(12) 



this gives a directly invertible relation for fi in terms of k 
and hence allows k to be estimated. 

For simplicity, we use the 'parabolic estimator' sug- 
gested by the cluster simulations of van Kampen (1998) 
which relates 7 to k via 



7 



II 



(13) 



As in T98, we adopt the value c — 0.7 corresponding to a 
profile which lies between that of a homogeneous sheet of 
matter (k = constant) and an isothermal profile {k cx , 
r the distance from the cluster centre). The magnification 
can therefore be written 



(14) 



where P = ±1 accounts for image parity on either side of 
the critical line implied by the modulus in equation (|l^. 
Rearranging this equation for k gives 



2c 

The second parity S is due to the parabolic nature of 
equation ( p^ which permits a second critical line and thus 
higher values of k. Since there is no evidence for the ex- 
istence of a double critical line in CL0024-I-1654, we will 
adopt 5 = +1 throughout this paper. 

Since the parabolic estimator requires only local k and 
7, it finds application to both radial and 2D magnification 
distributions in this paper. 



3.4. Self-consistent axi-symmetric mass estimator 

The second method which we use to estimate cluster mass 
is the so-called self-consistent axi-symmetric mass estima- 
tor introduced by T98. The non-local nature of this esti- 
mator allows the magnification equation ( |l^ ) to be solved 
for a self-consistent k and 7 radial profile. Although the es- 
timator is valid only for axi-symmetric mass distributions, 
it can be applied to data binned in any self-similar set of 
contours centred on the peak of the mass distribution. 
In a given annular bin i, the shear can be expressed as 



(16) 



where ki is the convergence in the bin and is the con- 
vergence averaged over the area interior to and including 
the bin. Substituting this into equation (|2|) gives for the 
magnification in bin i 



V^i;^ = (1-k,)(1-2k, 



(17) 



Dividing Ki into two terms, one for bin i and the other, 
which we denote 77^-1, for all interior bins, so that 



Vi-l 



i + l 



(18) 



allows rearrangement of equation ( p^ to give 
(* + l) 



S[i^-l-{^ + ^)^,_,f + A^V^lr^]'/^}. 



(19) 



The parities V and S have the same function as in Section 

Using equation (|l^), k is calculated iteratively. The 
only freedom is choice of tjq which as equations ( p^ and 
(|l8|) show, is 71, the shear in the first bin. To avoid non- 
physical solutions, 72 > V^i^ must be enforced leaving 
only a sensible range of values. As T98 discuss, the overall 
mass and shear profile obtained with the self-consistent 
axi-symmetric estimator proves to have only a minor sen- 
sitivity to the choice of 71 . 

3.5. Self-consistent 2D mass estimator 

The final reconstruction method we apply in this paper is 
that discussed by Dye & Taylor (1998). By pixellising the 
image into a rectangular grid of pixels, the components of 
the shear in any pixel n can be expressed as 



i- ((c2 + 1) - 5V(c2 + l)2-4c2(l-7'/i"i)) . (15) 7? 



1,2 



(20) 



where summation over index m is implied for all N pixels 
on the grid. The matrices Di and D2 are simple geometri- 
cal functions of the different combinations of positions of 
pixels ni and n (see Dye & Taylor 1998 for details). 

Writing equation (|l2|) in its pixellised form and substi- 
tuting the expression for shear in equation (^0|) gives the 
vector equation 



1 - 2k + kGk* - Vn" 







(21) 
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where V is again the image parity from Section 3.3, ^ is 



the iV-dimensional vector of pixelhsed inverse magnifica- 
tion values, is the transpose of the vector k of pixelhsed 
convergence values and 1 is the vector (1, 1, 1, ...). G is an 
N X N X N matrix with the elements 



(22) 



where 8ij is the Kronecker delta, and summation occurs 
only over indices p and q in equation (pT|). 

We solve equation (21) for k, using a hybrid Powell 
method provided by the NAG routine C05PCF. 



4. Cluster Model constraints from Number 
Counts 

Before application of the theory presented in Section ^ the 
properties of the background population must be under- 
stood. Our initial attention is turned toward the i? band 
population of background galaxies since this proves to be 
the most suitable for the application of lens magnification 
owing to its shallower number count slope. In this section, 
we use the radially binned number counts of the sample to 
fit an isothermal, power-law and NFW profile. We defer 
discussion of the U band sample properties until Section 



4.1. R band sample characteristics 

The characteristics of the R band background source 
galaxy population must be constrained before we can ap- 
ply our lensing analysis to the source counts. Specifically, 
the unlensed surface number density and the number 
count slope must be determined. 

T98 estimate the unlensed surface number density 
from deep number counts of field galaxies observed inde- 
pendently of their work. We opt for the alternative choice 
here of taking nu from our field where the lens effect of 
the cluster is expected to be small. This eliminates poten- 
tial biasing of the reconstructions which arise from differ- 
ing completeness characteristics between the lenscd field 
and the reference unlensed field (Gray et al. 2000). This 
is particularly important for the U band break search in 



Section 6.2. In an annular region centred on the cluster 
and bounded by the radial limits 120" < r < 180", we 
measure nn = 43 ± 2 arcmin"^ in the R. This value is 
assumed in our subsequent analysis. 

Ideally, this normalising annulus should be further 
away from the centre of the cluster than we have chosen. 
However, as we discuss in Section ||, the edge of our field is 
somewhat noisy due to the presence of some bright stars. 
This causes a drop in the number counts at larger radii 
so that normalising here would force a large over density 
of counts at medium radii and hence negative mass. An 
alternative explanation is that such an overdensity would 
be caused by contaminating cluster galaxies. In fact, it 
is most likely that both effects jointly contribute. Since 
our reconstructions include the uncertainty due to cluster 



contaminants, any bias as a result of normalising out this 
effect falls within our error budget. We therefore expect 
only a small bias to remain from the underestimation of 
background counts at large radii. 

As far as the number count slopes are concerned, the 
R band slope determined by Hogg et al. (1997) over 21 < 
i?AB < 25 is — 0.83, in agreement with Small et al. 
(1995) who measure [3b. — 0.80 over the same magnitude 
range. This is consistent with the slopes in the V and / 
from the Hubble Deep Field data of Pert et al. (1998) with 
f3v ^ 0.9 and /?/ ~ 0.8 over 23 < (F,/)ab < 26. However, 
over the fainter magnitude range 26 < {V,I)ab < 29, 
they measure a shallower slope of /3v,i < 0.5 and find the 
trend that flattening of the number count slope is more 
pronounced in shorter wavelength bands. The fact that 
this is seen at all in the / band therefore suggests that 
flattening most probably also occurs in R at these faint 
magnitudes. This is fortifled by the more recent R band 
counts of Metcalfe et al. (2001) from the William Herschel 
Deep Field, which suggest a slope of ~ 0.6 fainter than 
Rab ^ 26. 

Our R band observations tantalisingly extend to ap- 
proximately the depth whe re the apparent break occurs. 
We compare in Section O and |5.3| the difference between 
the results obtained using the single and dual slope mod- 
els. In the analysis hereafter, for the single slope model, we 
take = 0.80 and for the dual slope model, Pri = 0.80 
and f3ii2 — 0.5 with a break magnitude oi R — 26. 

4.2. Magnification profiles 

In the case of lens magniflcation, the least biased method 
of determining the best fit mass model is to fit depletion 
curves to the number count profile. In this way, the lensing 
signal is used in its purest form before potential biases are 
introduced by calculation of the k profile. 

To fit the depletion profile from a given mass model, 
its lens magnification must be determined. We choose to 
fit an isothermal sphere, a power-law mass model and a 
NFW profile, all of which have analytical forms for their 
magnification. In the case of the isothermal sphere, the 
mass profile is completely determined by its critical radius 

'^crit ; 



K{r) 



giving a magnification of 
fj,{r) = |1 - rcrit/r\~^ . 



(23) 



(24) 



The power-law model we choose is that of Schneider, 
Ehlers & Falco (1993) which gives a smooth, non-singular 
surface mass density distribution with a n profile given by 



k{x) — Kq 



2^2-4 ■ 



(25) 



Here kq is the peak surface mass density at the centre, ^ 
is the power-law slope and x — r/vc 
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radius. For this distribution to remain positive definite at 
large radii and so that it is a dechning function of radius, 
the vahd range of slopes isO<^<l. Atr:^ rcore, this 
physically corresponds to a range of mass models span- 
ning a homogeneous sheet = 1) through an isothermal 
sphere = 0.5) all the way up to k cx = 0). The 

magnification resulting from such a profile is (Schneider, 
Ehlers & Falco 1993) 



1 - 



Kq 



1 



(l + a;2)i 

Kq 



(l+a;2) 



,-{1 + (2^-1)^2} 



(26) 



Finally, the NFW model has a k profile given by 
(Bartelmann 1996) 



K{y) = 2ks— - 

- 1 

where 



(27) 



fiy) 



1 - 
1 





2 



tan 



tanh 




(28) 



and y — r/rs. The scale radius, r^, and surface density 
normalisation, k^, are free parameters of the model. The 
magnification of the NFW model comes from Jacobian 
determinant of the lens mapping for axisymmetric lenses 
(Schneider, Ehlers & Falco 1992); 



where 

m{y) = 2 



1 - 



m{y) 



dy 



m{y) 



(29) 



(30) 



4.3. Number count profile fits 

Using equation (^) for the single slope model (and neglect- 
ing the clustering term for now) with the three forms for 
/X from Section 4.2, we perform a fit to the observed de- 
pletion profiles. Figure ^ shows the R band radial number 
counts expressed as a fraction of the intrinsic background 
counts. The degree of obscuration by foreground objects 
is taken into consideration by adjusting in each bin. 
The la error bars plotted account for the uncertainty in 
Ufi and shot-noise. 

Superimposed on the observed R band number counts 
in Figure |^ are the isothermal, power-law and NFW pro- 
files obtained fro m fit ting to the first 9 data points. We 
discuss in Section 5A that the outer profile is affected by 
noisy features at the edge of the field of view. The final 
model fits therefore omit the last 3 points although we 
also fit to all 12 points to investigate the variation in fit- 
ted parameters. 




8 (arcmin) 

Fig. 7. Number of counts as a ratio of expected counts versus 
radial distance from cluster centre in tfie R band. Error bars 
account for the error in nu and shot noise. Superimposed are 
the isothermal (dashes, Vcrit = 15"), power-law (dot-dashes, 
Ko = 1.52, ^ = 0.61, reore = ll") and NFW (dots, r, = 1.75', 
Ks = 0.23) models fit to the first 9 points. 



For the isothermal profile, a best fit critical radius of 
rcrit — 15"±1" is obtained using either the first 9 points or 
all 12 points. Note that this is somewhat smaller than the 
estimate of Vcrit — 25" in the R band from ROl. This dis- 
crepancy arises partly from the different object selection 
and extraction criteria used by ROl, detailed in Section 
H, and also from our different radial bin width. The er- 
ror of ±1" we quote here is solely the error from the 
fit which does not include any uncertainty to allow for 
the choice of binning. Variation in bin width and also 
the radius attributed to a given bin adds further error. 
In Figure |^, the radial position of a given bin is taken as 
the radius which divides that bin into two equal areas. 
Since the isothermal fit is dominated by the radius of the 
first data point, adopting different binning strategies af- 
fects the fitted value of rcrit substantially. Measuring the 
variation in fitted values of rcrit with different binnings, 
we find that a further error of ±10" should be included to 
give an overall fit of rcrit — 15" ± 10". 

The observed radius of the large arcs, rare — 30", is 
also somewhat larger than our fitted critical radius. One 
explanation for this would be that the background popula- 
tion we select lies at a lower redshift on average than the 
lensed background galaxy forming the arc at z = 1.675 
(BOO). In fact, knowing the amount of mass contained 
within the arcs (see Section 5.3), one finds that a crit- 
ical radius of rcrit — 15" ± 10" corresponds to a back- 
ground population with a mean redshift lying within the 
range 0.40 < z^can < 0.82 (EdS or f7 = 0.3, A = 0.7). 
Note that although this range appears to permit only 
quite low mean redshifts, it depends sensitively on the 
confidence interval. For example, expanding the interval 
from 68% quoted above to merely 80% results in a range 
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0.40 < 

^mean ^ 1.19. Clearly, tlie large error on Tq^h pre- 
vents an accurate measurement of Zmean- 

However, considering the redshift survey of Cohen et 
al. (2000) in the Hubble Deep Field North region, this in- 
dicates a median redshift of Zmedian = O.TO^g'gJ in their 
deepest bin at -Rab — 24. Extrapolating their Zmedian ver- 
sus magnitude plot by a further 1.5 mags to coincide with 
our 3(7 detection limit and weighting by the expected num- 
ber counts over this interval, one obtains an approximate 
mean redshift expected for our background R sample of 
Zmean ~ 1-2. Therefore, although our lens-inferred mea- 
surement of Zmean IS a little low, it is not significantly 
inconsistent. 

The choice of bin width of 20" means that our abil- 
ity to constrain the small scale core radius in the power- 
law profile is very limited. In fitting to the power-law, we 
therefore hold Vcore = H" determined by the shear study 
of CL0024-I-1654 by TKD which has superior resolution 
in the centre of the cluster. Fitting the remaining two pa- 
rameters to the first 9 points then yields kq = 1-52 ± 0.20 
and ^ = 0.61 ± 0.11 with la errors accounting for the fit 
and binning variation. This compares to kq = 1.42 ± 0.21 
and ^ = 0.67 ±0.11 obtained when the remaining 3 points 
are included. 

Finally, fitting the NFW model to the first 9 points 
yields the parameters = 1.75' ± 1.02' and Kg — 0.23 ± 
0.08. Including all 12 points, these become = 2.50' ± 
1.32' and Ks — 0.19 ± 0.08. The large errors here reflect a 
strong degeneracy between Ks and r^. 

Although close to isothermal with a slope of ^ = 0.61, 
the power-law model provides a better fit to the number 
counts than the isothermal sphere. The NFW and power- 
law models are clearly very similar however the reduced 
of 0.41 ± 0.53 from the NFW model (assuming Gaussian 
statistics to obtain the error) quantifies the fact that it is a 
better fit to the data than the PL model with a reduced 
of 0.62±0.53 including the first nine bins. For comparison, 
the isothermal sphere model fits to the data with a reduced 

of 1.33 ± 0.50. This agrees with the findings of TKD 
who concluded that although the NFW profile predicts 
too much mass within the inner arc region, there is little 
to distinguish this from a power-law profile at larger radii. 
In fact, with all 12 bins in the fit, the PL model fares best 
with a reduced x^ of 1.08 ± 0.45 compared with that of 
the NFW model of 1.24 ± 0.45. 

Our magnification profile fits prove to be consistent 
with fits to shear-derived mass profiles from other work. 
TKD find that the mass profile of CL0024-I-1654 is best 
described by a power-law model with ^ = 0.57±0.02 and a 
central surface mass density of Eq — 7900± lOO/iM0pc~^. 
We can compare this central surface mass density with 
our estimate of by integrating our fitted power-law 
model for k over the disk enclosed by the observed arcs. 
Comparing this with the re al p rojected mass contained 
within this area from Section ^.3| allows the normalisation 
of the model to be calculated. This requires a central sur- 
face mass density of Eo = 9000 ± 1800/iMopc-2, slightly 



higher than the TKD estimate but consistent given the 
error budget. 

Similarly, BOO find that their azimuthally averaged 
mass profile for CL0024-I-1654 is close to the NFW predic- 
tion for massive clusters with an overdensity of 5c — 8000 
and Ts ~ 400/i~ikpc (NFW, Ghigna et al. 1998). Using 
the fact that 5c = i^s'^ck/ {ts Pcrit) (Bartelmann 1996) 
where EcR is the critical lens surface mass density (see 
for example Blandford & Narayan 1992) and pcrit is the 
critical density of the Universe, this converts to k,s ~ 0.2 
with a scale radius ~ 2'. The findings of this paper are 
therefore in good agreement with the NFW expectation. 

Finally, we have also applied the dual slope model to fit 
the three profile forms with magnification calculated using 
equation (|ll|). We find that this gives fitted model param- 
eters within a few percent of those presented above. As 
expected, this negligible difference indicates that the lens- 
ing signal in the R band is dominated by galaxies in the 
i? < 26 magnitude range where the steeper slope applies, 
despite fiux magnification of objects by CL0024-I-1654. 

5. Cluster Mass Reconstructions 

In this section, we present the results of the 2D and ra- 
dial mass reconstructions obtained using the R band back- 
ground galaxy sample. 

5.1. 2D mass maps 

Using the 12 x 12 grid of bins shown in Figure ^ and 
the associated R band obscuration mask, source numbers 
were binned across the field of view and used to calculate 
magnifications using equation (^. Applying both the local 
estimator and the self-consistent reconstruction method 
assuming single slope counts, we derive the mass maps 
shown in Figure ||. The top half of this figure shows the 
locally estimated n for comparison with the lower half 
showing the self-consistent mass distribution. 

Both maps exhibit a very similar mass structure 
demonstrating that the local estimator works well. The 
peak of the distribution in the self-consistent map is 
slightly higher than in the estimated map and is more 
distinct from the surrounding mass fiuctuations. This is 
especially true at the edge of the field where noisy fea- 
tures in the estimated map have been suppressed by the 
non-local solution provided by the self-consistent method. 
In particular, the features seen in the lower left and lower 
right corners due to a large degree of obscuration by the 
two brightest stars in the field (which is difficult to prop- 
erly account for) have been suppressed significantly. In 
addition, the large dip seen in the top left corner caused 
by an excessive number of background sources has been 
slightly reduced. 

The error map calculated from the width of the prob- 
ability distribution in equation (H) indicates that the sig- 
nificance of the peak in the self-consistent map is 6(T. If 
we include the extra uncertainty due to contamination by 
cluster objects at the predicted level of 19% as well as 
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Fig. 8. Mass reconstruction from the R band data. Top: 
Estimated k from parabolic estimator. Bottom: Iterated self- 
consistent K. Contours in both plots are separated by intervals 
of 5k = 0.1. North is up, East is left. 



K ~ 1.4. Furthermore, the local estimator agrees with the 
non-local estimator which would certainly not be the case 
if the bin resolution caused an effective smoothing of the 
profile. The cumulative mass measurements derived from 
the profiles are also therefore not affected. 

The topology of the structure seen in the mass maps 
is very similar to the distribution of cluster light (Figure 
^ and the cluster galaxy number density (Figure ||). In 
the mass map, we find a relatively compact core with ex- 
tensions toward the north-west and south-west (toward 
the top- right and bottom-right in Figure ^ respectively). 
In comparing these extensions with the light and num- 
ber density distributions, it is noticeable that the light 
traces the north-west extension more than the south-west, 
whereas the number density traces the south-west exten- 
sion more than the north-west. The explanation is that in 
the north-west, there are fewer galaxies, forming an almost 
distinct sub-clump, but they are relatively bright. In con- 
trast, the south-west extension harbours a higher number 
of fainter cluster members. 

Comparison of our mass map with the X-ray map of 
Soucail et al. (2000) again shows very similar structure. 
The X-ray contours show an elongation along the south- 
west direction while a separate peak of X-ray emission 
is seen toward the north-west, coinciding with the sub- 
clump of large cluster galaxies noted above. Soucail et 
al. suggest that the emission originates from one of two 
possible sources. The closest to the centre of the emission 
is a large cluster galaxy lying at a redshift of z = 0.4017 
while the other is a foreground star-forming galaxy ai z — 
0.2132. Given our detection of mass and the concentration 
of large cluster galaxies in this region, it seems plausible 
to explain at least part of the X-ray emission as being 
due to cluster gas. The fact that the X-ray map shows 
a distinct peak suggests that this is actually a separate 
mass clump rather than an extension. The shear map of 
TKD shows no evidence of there being such a sub-clump 
although BOO's map suggests that extra mass lies in this 
vicinity. Unfortunately, neither of these shear maps cover 
a wide enough field of view for this to be properly verified. 



the error in the unlensed surface number density, this sig- 
nificance drops to 4(7. Note that the peaks in both plots 
appear lower than the peak value implied by the radial 
analyses. This is a direct consequence of the relatively low 
resolution of the 12 x 12 grid of bins which has averaged 
over the central cluster region and effectively smoothed 
out the peak. This same reason describes the apparent 
discrepancy with the shear-derived mass maps of TKD 
and BOO, both of which cover a much smaller field of view 
(~ 1.5' X 1.5') and both of which imply a peak of k ~ 1.3. 

This resolution effect in the maps does not affect the 
shape of the radial mass profiles because it occurs on a 
scale smaller than the radial bin size. Section 5_^ verifies 
this since the profiles obtained agree with the power-law 
model which predicts a central surface mass density of 



5.2. Radial mass profiles 

In Figure ^ we show the radial mass profile calculated for 
the R band sources using both the local estimator and the 
axi-symmetric method assuming the single number count 
slope model. For calculation of the axi-symmetric profile, 
the val ue of the shear in the first bin is set to 71 = 0.25 (see 
Section 3^). The shaded regions in both plots show the la 
uncertainty accounting for shot noise, source clustering, 
uncertainty in background count normalisation and 19% 
cluster member contamination. Similar to the 2D mass 
distribution, we find that the last 3 data points in the 
axi-symmetric solution are slightly suppressed compared 
to the local estimator. 

The fitted profiles of Section O plotted in Figure || 
again show that the power-law and NFW models give a 
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Fig. 9. Top: Locally estimated radial mass profile for R band 
sources binned using annuli shown in Figure ^ The shaded 
region shows the la uncertainty which allows for shot noise, 
source clustering, uncertainty in background count normalisa- 
tion and 20% cluster contamination. Bottom: Axi-symmetric 
solution for the R band. Both plots show the isothermal 
(dashes), power-law (dot-dashes) and NF W ( dots) models fit- 
ted to the first 9 data points (see Section 4.3). 



better fit to the results. Both profiles exhibit an excess 
of mass at 3" < r < 4". This is attributed to the noisy 
features seen in the 2D mass plots at the edge of the field, 
and just pushes the measured mass at this radius to a 
value inconsistent with the NFW and power-law fit. 

5.3. Cumulative mass profile 

To convert from k to real projected mass, we normalise to 
the amount of mass contained in the disk traced by the 
observed arcs. Since the enclosed mass can be calculated to 
a high accuracy, normalising in this way provides a much 
more reliable scaling than using a mean redshift estimated 
for the background population. Using BOO's measurement 
of the redshift of the arcs at z = 1.657, the mass within 
a radius of rare = 30" from the centre of the cluster can 
be calculated to be 1.19 x IO^'^H'^Mq for an Einstein de 
Sitter (EdS) Universe. The cosmological dependence of 
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Fig. 10. Cumulative projected mass calculated by normalising 
to the amount of mass enclosed within the circle traced by the 
observed arcs. Also plotted is the NFW model of Section 4.3 
(dashes) and the B band results of van Kampen (1998) (dot- 
dashes) . 



this result is weak to the extent that for an = 0.3, A = 
0.7 cosmology, this enclosed mass increases by only 5%. 
This scaling applies generally hence all aperture masses 
quoted hereafter assume an EdS cosmology. 

Scaling our radial k profile with this normalisation, we 
sum the mass in each bin to produce the cumulative pro- 
jected mass profile in Figure |l^. Out to a radius of 2.9' 
(0.54 Mpc for EdS), we measure a projected mass of 
(8.1±3.2) X 10^^/i~^Mq. The la error here again includes 
shot noise, source clustering, uncertainty in background 
count normalisation and contamination from cluster mem- 
bers at the 19% level. If we apply the dual slope model, 
we find a mass of (8.2 ± 3.2) x IO^^Ii'^Mq, re-iterating 



the statement made in Section 4.3 that the lensing signal 
in the R band is dominated by galaxies in the i? < 26 
magnitude range where the slope Pr — 0.80 applies. 

Approximating the integrated NFW profile, the pro- 
jected mass within a radius R scales as M(< R) = 
2.9 X 10^'^{R/l')^-^'"-^^s{R/i')h-^MQ. With this scaling 
relation for the cumulative mass, we can readily com- 
pare our results with other authors. Considering exist- 
ing lensing measurements for the moment, the measure- 
ment by BOO of the mass contained within the arcs ob- 
viously agrees with our estimation by default since this 
is essentially the mass we normalise to. Using the red- 
shift of the arc from BOO, TKD's estimate of the clus- 
ter mass contained within a radius of 107ft.^^kpc is 
(1.553 ± 0.002) X IQI^/i-IMq. At this radius, our NFW 
estimate for the enclosed mass is (1.3 ±0.3) x lO^^/i~^M0. 
Extending to slightly larger radii, the magnification anal- 
ysis of van Kampen (1998) finds a projected mass of 
(6.7 ± 2.0) X IO^^'/i-^Mq within SOO/i'^kpc compared 
to our NFW mass of (5.0 ± 1.7) x lO'^'^h-^MQ. On even 
larger scales, the weak shear study of Bonnet, Mellier& 
Fort (1994) measures a mass of (1.6 ± 0.4) x IO^^/i-^Mq 
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within R = 1.5 h ^ Mpc compared with our prediction of 
(1.7 ±0.7) X IO^'^H-^Mq. 

In terms of mass-to-Hght, using the surface flux den- 
sity of cluster galaxies plotted in Figure ^, we find a ratio 
of M/Ls = 330 ± 30 (total mass/total light in B) within 
a radius of 0.5'. This error accounts for uncertainty in the 
photometric zero point but does not reflect the fact that 
this is effectively an upper limit as a result of potentially 
missing cluster galaxies in our selection process (Section 
2.3). Again, since our mass within this radius agrees with 
that of BOO, we find an almost identical result to BOO's 
result of M/Lb = 320 ± 30. On a slightly larger scale, we 
measure M/Lb — 480 ± 180 within 1' where now the er- 
ror includes additional uncertainty from the mass error. 
Finally, out to 2.9', we find M/Lb = 470 ± 190. These lat- 
ter two measurements are consistent with a mass-to-light 
ratio which increases on larger scales although a scale- 
independent ratio, as found by Kochanski, Dell' Antonio 
& Tyson (1996), remains plausible given the errors. 

In Figure |l^, we also plot the only other magnifica- 
tion based cumulative mass profile of CL0024-I-1654 which 
exists to date; that of van Kampen (1998). This was de- 
rived from the B band observations of the cluster by Fort, 
Mellier & Dantel-Fort (1997). We choose the B band data 
in favour of the / which suffered contamination from an 
unusual speckle beam pattern coming from stray light 
in the telescope optics (Fort 2001, private communica- 
tion). This resulted in an over-exaggerated zone of deple- 
tion in the / band radial number counts. Notice that van 
Kampen's B band mass profile agrees very well with our 
data points in terms of its shape but that its normalisation 
is higher by ~ 50%. 

Turning to alternative methods for the determination 
of cluster mass, Czoske et al. (1999) using galaxy dynamics 
from 227 spectra measure the 3D mass of CL0024+1654 to 
be (1.4±0.3) X lO"/i-^M0 within a radius of 500 h-^ kpc. 
Converting this to a projected mass using the isother- 
mal sphere model they assume, yields a value of (2.2 ± 
0.5) X 10^'^h~^MQ. This agrees with the X-ray tempera- 
ture measurements by Soucail et al. (2000) who measure 
a projected mass of 2.3tli X IO^^'/i-^Mq within the same 
radius. Our scaling relation gives a projected mass en- 
closed by this radius of (8.0 ±3.1) x IO^'^H-^Mq. It should 
be noted that the later redshift measurements of Czoske 
et al. (2001a) indicate that CL0024+1654 is a perturbed 
system. This implies that mass measurements of the clus- 
ter which assume dynamical equilibrium are probably not 
reliable. We discuss the significance of this in Section |^. 



6. The U band galaxy population 

The role of the U band population thus far has been 
to assist in the selection of R band background galaxies 
through colour cuts. In this section, we test the feasibility 
of using the U band for detection of lens magnification 
through number counts. We wish to re-examine the claim 
made by ROl that the number count slope of the back- 



ground field galaxy sample becomes flatter than the lens 
invariant slope, (3 = 1, at faint magnitudes. 

6.1. The U band number count slope 

Current evidence regarding the notion of a break in the U 
band field galaxy number count slope at faint magnitudes 
is tenuous. At magnitudes brighter than Uab — 25.5, 
there is good convergence. All studies agree that at these 
brighter magnitudes, the number count slope is relatively 
steep. Pozzetti et al. (1998), in analysing the Hubble Deep 
Field North (HDF-North), find that for 23 < Uab < 25.5, 
(3u = 1.00 whereas Hogg et al. (1997) conclude that 
(3if — 1.17 for Uab < 25.5 from ground based observa- 
tions. 

At fainter magnitudes, Pozzetti et al. (1998) measure 
a much fiatter slope of Pu = 0.34 at Uab = 25.8 from the 
HDF-North. Similarly, at the same magnitude, Metcalfe 
et al. (2001), using a combination of the groimd based 
William Herschel Deep Field and the HDF-North and 
HDF-South, report a slope of f3u = 0.38. This appears 
to contrast with the results of Volonteri et al. (2000) who 
measure a slope of = 1 up to Uab — 26.5 with the 
HDF-South. 

The depth of our observations in the U band extend 
to approximately the suggested U band break magnitude, 
which at first sight, appears to extinguish any hopes of 
pushing the search beyond current limits. However, we are 
in a fortuitous position for two reasons. The first is that 
our analysis merely rests upon the relative depletion of 
galaxy numbers. As Gray et al. (2000) demonstrate, this 
means that provided the completeness characteristics of 
the unlensed reference number counts are the same as the 
lensed field, the completeness function drops out of equa- 
tion d^). Our reference counts are taken from the edges 
of our cluster field so by definition have the same com- 
pleteness characteristics. In the same way that the com- 
pleteness function drops out of equation (^, it also drops 
out of equations ( 11 ) and (^) . This is an important result 
and means that our break search is not affected by incom- 
pleteness. We therefore need not concern ourselves with 
completeness corrections such as those applied to existing 
traditional number count measurements which may cause 
a potentially large uncertainty at faint magnitudes. 

Our second fortuity lies in the fact that our obser- 
vations are of a cluster-lensed field, enabling us to exploit 
the effect of flux magnification on our background sample. 
Taking the NFW model fit from Section 4.3, this predicts 
that within a disk of diameter 0.8' centred on the clus- 
ter, lens magnification by CL0024+1654 pushes the effec- 
tive limiting magnitude of our observations by Am > 2.5. 
Within a disk of diameter 2', this becomes Am > 0.8. 
We are therefore able to use our dataset to search for a 
break in the U band number counts some way beyond the 
physical limiting magnitude of our data. 

Note that this does not contradict the discussion ear- 
lier that allowing for a dual slope in the R band lensing 
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analysis causes little difference to the results obtained. It 
is true that both the U and R band samples extend to 
approximately where the break magnitude is thought to 
lie and it is also true that lens magnification acts on both 
samples equally to increase the effective depth of each. 
In the case of the R band sample however, the determi- 
nation of mass depends on the joint contribution from 
both the steep and shallow number count slopes. In the 
search for the U band break as we discuss in detail be- 
low, these contributions are considered separately so that 
the more numerous galaxies in the steep slope region do 
not dominate those beyond the break magnitude. In ad- 
dition to this effect, the expected U band change of slope 
of /3(7 ~ 1 ^ 0.3 is more severe than the expected R band 
change of jSu ~ 0.8 — > 0.5. This makes detection of any 
faint slope depletion easier and hence further increases the 
confidence limits we can place on our search results. 

6.2. Testing for a U band break 

We use two number count models to search for a break 
magnitude. The first is a modification of the single slope 
number count scenario to allow the slope to smoothly flat- 
ten off beyond the break ma gnit ude. The second uses the 
dual slope model of Section |3.2| , allowing the faint slope 
to vary. 

6.2.1. Model 1 

The first model generalises the single slope number count 
model to include an arbitrary cut-off at Uq; 



I ^ IqO.40uAU 



1 + /iA-'100-4/3l'AC/ 



(31) 



where AU = Uum — Uo and Uum is the limiting magnitude 
of the observations, Uum = 25.7. Well below the break 
scale, Uq <C Uum, or for weak lensing « 1, this reduces 
to the usual scaling for lensing of a single slope number 
count distribution in equation (^. Above the break scale, 
equation ( |3l| ) tends toward the scaling n' = n/i^^, for a 
completely fiat number count distribution. 

We take the best fit NFW profile determined from the 
R band data to provide the magnification in equation (|3l|). 
The free parameters of the model are therefore the break 
scale, Uq, and Pu which we fit by maximising the following 
likelihood function: 



n 



exp 



nu,i 




(32) 



Here, n'jj ^ is the number of U band objects observed in bin 
i, nu,i is the number expected in the absence of lensing 
and pu,i is the number predicted by the NFW profile fit- 
ted to the R band data. The product acts over all annuli. 
Uncertainty due to cluster and foreground galaxy contam- 
ination is incorporated through the quantity q. The fore- 
ground galaxy component of Ci is set at 3% of nu^i for 
each bin, in accordance with the contamination fraction 





Fig. 11. Top: x" 
slope Pu and the break scale Uo assuming the number counts 
flatten off completely at faint magnitudes. Bottom: distri- 
bution of the faint number count slope f3u2 and the break scale 
holding the bright slope at Pui ~ 1.07. AU contours are sepa- 
rated by Ax^ ~ 1. 



discussed in Section 2.4.1 . The cluster component of 
varies with choice of CLF. We therefore average over dif- 
ferent realisations of the CLF, weighting by the probabil- 
ity distribution obtained from fitting to our known U band 
cluster counts (see Section 2.4.2| ) . As in Section the ra- 
dial variation of cluster number count density is taken to 
be k/r where k is set by normalising to the contamination 
determined for each CLF realisation. 

The top half of Figure |ll| shows the contours ob- 
tained from the likelihood distribution of equation (p2|). 
With a 95% confidence, our data rules out a complete 
fiattening of the U band counts with a break magnitude 
brighter than Uab = 27.3. Furthermore, we measure the 
U band slope as (3u — 1.07±0.06. This is very close to the 
lensing invariant slope P — 1 which causes the increase in 
surface number density due to flux magnification to ex- 
actly cancel the dilution caused by magnification of their 
inter-spacing. The claim made in Section ^ that the U 
band population is unsuitable for the detection of magni- 
fication from number counts is thus substantiated. 
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Fig. 12. Radial number count profile in the U band. 
Superimposed are the best fit isothermal (dashed line) , power- 
law (dot-dashed line) and NFW (dotted line) models using the 
number count slope Pu ~ 1.07. Error bars account for the error 
in 71(7 and shot noise. 



Figure shows the U band radial number count pro- 
file observed. In this plot, the model fits determined in 
Section 4.3 are used to predict counts assuming jSu 



1.07. 

The predictions and measurements are clearly consistent 
with each other. There is little distinction between any 
of the predicted profiles due to the near lensing invariant 
slope. This serves to demonstrate the insignificant deple- 
tion signal imposed on the number counts in the U band. 

6.2.2. Model 2 

The second method we use to search for a break assumes 
the dual number count slope model described in Section 
3.2. For the purpose of our analysis in this section, we 
hold the bright slope at the value determined previously, 
Pui = 1.07. The faint slope (3i/2 and the break magnitude 
are allowed to vary in our minimisation. 

The lower half of Figure |ll] shows the results of this 
minimisation. Slightly more relaxed than the previous re- 
sult, the data rule out a complete flattening {Pu2 = 0) 
brighter than Uab = 26.6 at a confidence level of 95% . 
As the faint slope is allowed to steepen, this limit slips 
to brighter magnitudes as a result of the degeneracy be- 
tween both parameters (Pozzetti at al. 1998). However, 
a slope of f}u2 = 0.4, the value reported by Pozzetti et 
al. (1998) and Metcalfe et al. (2001) as being applica- 
ble fainter than Uab — 25.5 — + 26, can be ruled out to 
Uab < 26.4 with 95% confidence. Pushing this limit fur- 
ther still, the contours in Figure |l^ show that a slope even 
as steep as /3y2 = 0.8 can be ruled out with 95% confi- 
dence to Uab < 26.0. Finally, at (3u2 — Pui, where the 
degeneracy between the faint and steep slope is 100%, the 
break magnitude cannot be constrained at all. 

We note that these results are also inconsistent with 
the findings of ROl who suggest a combined slope (ie. the 
effective single slope causing the same level of depletion 



as a dual-slope model, steeper than the faint end slope) of 
13 ~ 0.5. 

7. Summary and Discussion 

From U and R band observations of the cluster 
CL0024-I-1654, we have selected a background sample of 
galaxies in both bands using colour information. The clus- 
ter member contaminants identified by matching with the 
Czoske et al. (2001a) redshift survey of the field have been 
discarded from the U band. We have computed the clus- 
ter luminosity function of CL0024+1654 in both bands. 
This has been used in conjunction with the U and R field 
galaxy luminosity function from the CN0C2 survey (Lin 
et al. 1999) to estimate foreground and cluster contami- 
nation of our background samples. In the R band, we esti- 
mate a total contamination of 21% in contrast to a smaller 
12% in the U, where the cluster component is calculated 
from our own fitted CLF. 

The shallower number count slope observed in the R 
band sample makes this suitable for an investigation into 
the lens magnification induced depletion of background 
number counts. Depletion in this sample has indeed been 
detected and used to measure the radial and two dimen- 
sional distribution of mass in the cluster. Out to a ra- 
dius of 0.54 /i^^ Mpc, we measure a total projected mass 
of (8.1 ± 3.2) X IO^'^/i-^Mq where the la error includes 
shot noise, source clustering, uncertainty in background 
count normalisation and contamination from cluster and 
foreground galaxies. We find that this result increases by 
merely 2% when we allow for a flattening in the R 
band number count slope at Rab = 26. Claims of such 
a change in slope have been reported by several authors 
(eg. Pozzetti et al. 1998; Metcalfe et al. 2001), however 
allowing for it in our depletion analysis makes little differ- 
ence due to the lensing signal being dominated by galaxies 
brighter than the break. 

Converting the R band flux from selected cluster galax- 
ies to the B band, we find a mass-to-light ratio of M/L^ = 
330 ± 30 inside an aperture of radius 0.5' centred on the 
cluster. Since our selection process will have inevitably 
missed cluster galaxies, neglecting their contribution to 
the total luminosity, this is effectively an upper limit. On a 
slightly larger scale, we measure M/Lb — 480±180 within 
1' and also M/Ls = 470 ± 190 within 2.9' (0.54 /i"! Mpc). 
These latter two measurements are consistent with a mass- 
to-light ratio which increases on larger scales although a 
scale- independent ratio remains plausible given the errors. 

We have compared the observed radial depletion with 
that expected from an isothermal, power-law and NFW 
mass profile by fitting the predicted magnification of 
each. The NFW model provides the best fit to our data 
however there is little distinction between this and the 
power-law model. This agrees with the findings of Tyson, 
Kochanski & Dell' Antonio (1998) as does the mass 
of 1.3 X 10^'^h~^MQ we measure within the disk de- 
scribed by the large arcs compared to their estimate of 
(1.553 ± 0.002) X IO^/i-^Mq. Approximating the fitted 
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cumulative NFW model, the projected mass contained 
within a radius R scales approximately as M(< R) — 
2.9 X IQ^^iR/Vf-^-^-^'^^^^/^'^h-^MQ. Together with the 
cluster galaxy light and number density distribution, our 
2D mass maps suggest the existence of a separate sub- 
clump of mass just north-west of the cluster centre. This 
claim is strengthened by the detection of distinct X-ray 
emission in this area (Soucail et al. 2000) as well as the 
suggestion of extra mass from the shear mass map of 
Broadhurst et al. (2000). 

In fitting the number counts to the isothermal sphere 
mass model, we find that the fitted critical radius is 
smaller (although not significantly given the error) than 
the radius of the circle traced by the observed large arcs. 
This might be due to the fact that our selected back- 
ground population of sources lies at a lower mean red- 
shift than the lensed galaxy forming the arc at z = 1.675 
(Broadhurst et al. 2000). The mean background source 
redshift in an 80% confidence interval is inferred to be 
within 0.40 < Zmean < 1.19. The depletion analysis of 
CL00244-1654 by Fort, Mellier & Dantel-Fort (1997) in 
the B band finds a relatively wide depletion profile (see 
Section 5.3 regarding contamination of their / band data). 
Simultaneously fitting to a maximum of five isothermal 
sphere models, they find that a combination of various 
background galaxy populations is required to explain the 
wide zone of depletion. In B, this results in the mean 
background source redshifts ranging over z = 0.9ta l to 
z = S.OtoJ with 42% of sources lying at z = 0.9. This 
is a little surprising given that we find an adequate fit to 
our observed depletion profile assuming the existence of 
only one critical line. The discrepancy would perhaps be 
at least partially resolved by fitting the Fort, Mellier & 
Dantel-Fort data to an NFW profile since this accommo- 
dates a larger zone of depletion compared to the isother- 
mal sphere model. 

Comparison of our results or indeed any of the existing 
lensing results with those measured using cluster galaxy 
dynamics or X-ray temperatures shows a large discrep- 
ancy. We predict approximately 3.5 times as much pro- 
jected mass as the X-ray and dynamical measurements 
imply. Soucail et al. (2000) discuss that one explanation 
for this discrepancy may come from the inability to cor- 
rectly measure lens shear combined with a lack of accu- 
rate knowledge of the background source redshift distri- 
bution. While this bears some truth in general, it is not 
the case here. The determination of the arc redshift by 
Broadhurst et al. (2000) makes the measurement of pro- 
jected mass contained within the arc radius a robust one; 
this is a result which does not depend on the detection of 
shear through weak lensing and also knowing the redshift 
of the arced galaxy breaks the source redshift degener- 
acy. The projected mass contained within the arc radius 
according to the X-ray measurements of Soucail et al. is 
(0.5 ± 0.3) X lO^^h-^MQ. This is still a factor of nearly 
three times smaller than the arc predicted mass. 



An alternative scenario which provides a more satisfy- 
ing answer to the evidence gained thus far is alignment of 
additional foreground and/or background mass along the 
line of sight to CL0024+1654. The Czoske et al. (2001a) 
redshift survey identifies a group of galaxies lying just in 
front of the cluster as well as a pair of groups lying be- 
hind it. During the course of preparing this paper, new ev- 
idence emerged regarding the perturbed galaxy dynamics 
of CL0024-I-1654. Comparing numerical simulations with 
their redshift measurements, Czoske et al. (2001b) con- 
clude that the cluster could in fact be the result of a high 
speed collision of two smaller clusters along the line of 
sight. This would certainly explain why lensing observes 
a factor of three times as much mass. It also explains why 
early measurements of velocity dispersion were very large 
(see Dressier & Gunn 1992 and references therein); these 
were based on ~ 10 times fewer redshifts and hence failed 
to properly resolve the cluster's true dynamical state. 

Finally, the U band selected background sample indi- 
cates a near-lensing- invariant slope oi f3 — 1.07 ±0.06 and 
hence does not exhibit any noticeable sign of depletion. We 
have used the fact that lens magnification of our U band 
sample allows us to search deeper than its physical limit- 
ing magnitudes thus facilitating the search for a change of 
slope in the U band number counts at faint magnitudes. 
Up to Uab < 26.6, we can rule out a complete flattening 
to 95% confidence. Furthermore, we can rule out the ex- 
istence of a change of slope of /3 = 1 ^ 0.4 reported by 
Pozzetti et al. (1998) and Metcalfe et al. (2001) with a 
confidence of 95% up to Uab < 26.4. 

Our findings also contradict the results of ROl who 
claim to have measured U band depletion and hence slope 
flattening. Although we started with the same observa- 
tions as ROl, we applied more stringent criteria to the ob- 
ject extraction and carried out a thorough U band break 
analysis which allowed for contamination by cluster and 
foreground galaxies. It is intriguing that the HDF- South 
U band number counts reported by Volonteri et al. (2000) 
also support our result. They too find no evidence of flat- 
tening in the U band counts up to UABsimeq26.5. This will 
most likely remain an unresolved issue until deeper obser- 
vations, particularly in a cluster environment where lens 
magnification provides natural assistance, are obtained. 
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